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Abstract 

We study the properties of light mesons in the scalar, pseudo-scalar, and vector channels within 
the light-front quantization, by using the (one flavor) Nambu-Jona-Lasinio model with vector 
interaction. After taking into account the effects of chiral symmetry breaking, we derive the bound- 
state equation in each channel in the large N limit {N is the number of colors) , which means that 
we consider the lowest qq Fock state with the constituent quark and antiquark. By solving the 
bound-state equation, we simultaneously obtain a mass and a light-cone (LC) wavefunction of the 
meson. While we reproduce the previous results for the scalar and pseudo-scalar mesons, we find 
that, for a vector meson, the bound-state equations for the transverse and longitudinal polarizations 
look different from each other. However, eventually after imposing a cutoff which is invariant under 
the parity and boost transformations, one finds these two are identical, giving the same mass and 
the same (spin-independent) LC wavefunction. When the vector interaction becomes larger than 
a critical value, the vector state forms a bound state, whose mass decreases as the interaction 
becomes stronger. While the LC wavefunction of the pseudo-scalar meson is broadly distributed in 
longitudinal momcntmn (x) space, that of the vector meson is squeezed around x = 1/2. 
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1 Introduction 



Understanding the structure of light mesons at low energies is one of the most difficult problems in 
QCD, which requires nonperturbative studies about three entangled issues: chiral symmetry breaking, 
confinement of quarks and gluons, and relativistic bound-state physics. The light-front (LF) quanti- 
zation is now widely accepted as the most transparent formalism for the third issue, the relativistic 
bound-state physics [1], and thus, if the other two are also accommodated by the LF formalism, then 
it can be a very powerful tool for the problem. This expectation has been shared by many people 
with continuous interests and efforts in the LF formalism. The present paper is one of such attempts 
towards establishing the methods for describing the chiral symmetry breaking on the LF. 

The advantage of the LF formalism in treating the relativistic bound-state physics is closely related 
to the fact that the vacuum in this formalism is very simple, and is actually the Fock vacuum itself 
even for an interacting system. Roughly speaking, this property allows us to construct a Hamiltonian 
equation i^LC \oi) = E \a) with only a few Fock states just like in the Tamm-Dancoff approximation. 
By solving this equation, one can simultaneously find the wavefunction of each Fock state and the 
eigen-value (mass) of the eigen state. This is certainly the unique merit of the LF formalism. On the 
other hand, the other two properties, chiral symmetry breaking and confinement, are usually believed 
to be related to the nontrivial vacuum structure. This does not mean that the LF quantization is 
not useful for the problem mentioned above, but rather means that one has to discover new ways of 
formulating the physics of confinement and chiral symmetry breaking within the context of the LF 
quantization. This motivated several efforts towards setting up the problems in QCD, and actually 
many new interesting ideas and techniques were generated in such activities (see for example, Ref. [2]). 
However, we have not reached at a satisfactory description for these problems. 

Meanwhile, there has been considerable progress in describing the spontaneous symmetry break- 
ing in scalar theories within the LF quantization [1]. It was recognized that the longitudinal zero 
mode of a scalar field plays a crucial role. The scalar zero mode is not a dynamical mode in the LF 
quantization, and subject to a non-linear constraint equation, called the zero-mode constraint. Spon- 
taneous symmetry breaking can be achieved through the nonperturbative solution to the zero-mode 
constraint, which brings in a nonzero vacuum expectation value to the scalar field. Moreover, based 
on the analogies with such achievements, even the chiral symmetry breaking in simple fermionic mod- 
els (without gauge fields) were described within the LF quantization [3, 4] (See Ref. [5] for a review 
about these approaches based on the canonical LF quantization. For an alternative treatment based 
on the path-integral approach, see Ref. [6]). Therefore, the region governed by the LF formalism is 
now slowly expanding, while the chiral symmetry breaking in real QCD is still not under control (see, 
however, Ref. [7]). 

Indeed, within the Nambu-Jona-Lasinio (NJL) model, which is sometimes considered as an effec- 
tive theory of QCD [8, 9], it has been found that, still with the trivial Fock vacuum, chiral symmetry 
breaking is described in such a way that one selects an appropriate Hamiltonian depending on the 
phases of the symmetry [4] . In the NJL model, different Hamiltonians are originated from different 
solutions to a constraint equation, which exists only in the LF formalism. The "bad" component of 
the spinor is not a dynamical variable and is subject to a constraint equation like the zero mode in 
scalar theories. This "fermionic constraint" is nonlinear in the NJL model and leads to the "gap equa- 
tion" for the chiral condensate. When the coupling constant is larger than the critical value, the gap 
equation has a non-zero solution even in the chiral limit. This means that the fermionic constraint 
allows for "symmetric" and "broken" solutions corresponding to those of the gap equation. If one 
selects the "broken" solution, and substituting it to the canonical Hamiltonian, one obtains the "bro- 
ken" Hamiltonian. This governs the dynamics in the broken phase and is completely different from 
the Hamiltonian with the " symmetric" solution. In Ref. [4] , two of us solved the fermionic constraint 
by using the 1/N expansion, and obtained the Hamiltonians in both symmetric and broken phases. 
They also solved the bound-state equations for the scalar and pseudo-scalar mesons and obtained 



3 



their light-cone (LC) wavefunctions and masses, as well as the PCAC and Gell-Mann, Oakes, Renner 
(GOR) relations. 

The main objective of the present paper is to extend the analysis of Ref. [4] to the vector channel 

by using a similar model, and investigate the properties of a vector meson. A part of the results was 
already reported without detailed derivations [10]. Thus, in this paper, we will re-derive all the results 
in a self-contained way as much as possible. We will also present a slightly deeper understanding of 
the properties of scalar and pseudo-scalar mesons. 

It will be helpful to briefly summarize the main results of the present paper. After lengthy 
calculations, we are able to derive qq bound-state equations in the pseudo-scalar (vr), scalar (cr), and 
vector (V) channels. By solving these equations, we will eventually obtain the LC wavefunctions and 
masses of the mesons. Remarkably, it turns out that, due to the contact interaction of the NJL model, 
the spin-independent part of the LC wavefunction for each meson {i = tt, cr, V) has a very simple form 
(see Eq. (5.1)): 

where x is the fraction of longitudinal momentum carried by a quark, k± is the relative transverse 
momentum between the quark and the antiquark, M is the constituent quark mass, and rui is the 
mass of the meson i. This result is common for all the mesons discussed in the paper. The difference 
among the LC wavefunctions is visible only through the value of the meson mass m,. For example, in 
the chiral limit, we will find a relation = m-,^ < niy < nia = 2M. Namely, the LC wavefunction of 
the pseudo-scalar has a broad distribution in x space, while that of the scalar meson is rather peaked 
around the mean value x = 1/2, and the vector meson comes in between these two (see Fig. 3). 

The paper is organized as follows. In the next section, we will define the NJL model with the 
vector interaction. The vector interaction is added to have a bound state in the vector channel. We 
also show the explicit form of the fermionic constraint of this model. In Sect. 3, we solve the fermionic 
constraint by using the 1/A^ expansion, and derive the LC Hamiltonian to the first nontrivial order. 
This allows us to derive the bound-state equations of light mesons, which is discussed in Sect. 4. To 
the leading nontrivial order of the 1/N expansion, mesons are described by a quark-antiquark state, 
and the bound-state equations can be restricted to the lowest qq Fock states. In Sect. 5, we solve the 
bound-state equations in scalar, pseudo-scalar, and vector channels, and obtain the masses of mesons 
and the LC wavefunctions as already discussed above. We will show that, when we treat the bound- 
state equations, we have to be extremely careful about the regularization scheme. Otherwise, the 
equations for transverse and longitudinal polarizations of the vector meson look differently. Summary 
and discussions are given in Sect. 6. Some detailed calculations, as well as the notation we use are 
presented in Appendix. 

2 The model 

In this section, we define the NJL model we use, and discuss its particular properties when it is 
quantized on the light front. We also clarify our strategy which we take in the present paper. 

2.1 Nambu— Jona-Lasinio model with vector interaction 

The simplest NJL model is given by 

£S = ^ [i^a^af + i^an5^af] , (2-1) 

where is a spinor in the fundamental representation of "color" SU(A'') group (a is a color index: 
a = 1,---,N). The scalar interaction Cs is constructed as "color" singlet. In what follows, we 
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suppress the color index a for notational simplicity, which however does not cause any confusion 
because, throughout the paper except this section, we always treat color singlet objects (i.e., bilocal 
operators to be defined below). We consider the case with one flavor for simplicity. Generalization to 
the case with multi flavors, which is necessary for phenomenological applications, is more complicated 
but is straightforward. Recently, this model was analyzed in the light-front quantization by two of 
the authors [4] (see also Ref. [3]) and they computed masses and light-cone wavefunctions of scalar 
and pseudo-scalar mesons to the first nontrivial order of the 1/N expansion. One can of course apply 
the same procedure for vector states, but we know that the simplest NJL model defined by £njl does 
not allow for a bound state in the vector channel [8, 11]. Vector states start to bind if one adds the 
vector interaction so that the attractive force between a quark and an antiquark in the vector channel 
becomes stronger. Therefore, we include the vector interaction minimally by adding the following 
interaction: 

Notice that this interaction does not change the nontrivial "vacuum structure" (more precisely, the 
value of chiral condensate) which is caused by Cs- This is because wc have taken the normal order in 
(2.2) which is defined with respect to the Fourier modes of the fermion field (see below for the details). 
On the other hand, the bound-state equations for scalar and pseudo-scalar mesons will be modified 
due to the vector interaction. 

It would be very helpful to summarize here the issues to be discussed in the present paper and 
our strategy against them. There are mainly three issues: (i) chiral symmetry breaking, (ii) qq bound 
states in the scalar and pseudo-scalar channels, and (iii) qq bound states in the vector channel. Below 
we briefly explain the roles played by the interactions jCs and jCy for these problems. The precise 
meaning of all the procedures and approximations will be clarified in the forthcoming sections. 

(i) Chiral symmetry breaking is generated by the interaction jCq, and is not affected by the 
inclusion of the vector interaction jCy (by construction) to the leading order of the 1/N expansion. 

(ii) Bound states in the scalar and pseudo-scalar channels are possible^ in the presence of 
the scalar interaction Cs- Thus, in this paper, we will study them with the Lagrangian (2.1). 
Although the vector interaction Cy will modify the scalar and pseudo-scalar bound states, we 
shall ignore the effects of £v expecting that the scalar interaction Cs gives the dominant effect. 

(iii) Bound states in the vector channel are not formed with Cs alone, but can exist with the 
help of Cy. Thus, to study the bound-state physics in the vector channel, we need to include the 
effects of Cy. In the actual calculation, we will take the leading contribution of the eigen-value 
equation with respect to Gy which is enough to form a bound state. 

Therefore, in the following analysis, inclusion of the vector interaction (2.2) is only relevant for the 
bound-state physics of the vector state. 



2.2 A constraint equation for the bad spinor component 

One of the main sources of complexities in the light-front quantization is the presence of constraint 
equations which appear only when we treat the light-cone variable x"*" as the evolution time. In a 
system involving a fermion field, one always has a constraint equation for the "bad" component of 
the spinor ^_ (where V'± = = 57^7=*=* and V'-i- is refereed to as "good", see Appendix A 

for our notations) because the kinetic term i^i^"^ is expressed as iip^^d+ip^ + i^/^lc?-^- and 
is now a spatial derivative d- = d/dx~. This constraint equation which we call the "fermionic 

"^Precisely, to the first nontrivial leading order of the 1/N expansion, the scalar meson appears as a bound state only 
in the chiral limit. Besides, even with the vector interaction included, this result is the same because it affects only the 
pseudo-scalar channel. 



5 



constraint", is easily obtained if we multiply the Euler-Lagrange equation by the projector A_|_ from 
the left {A+j-d-'i> = -f-A.d-'if = -f-d-ip-): 



= 0. 



Without the vector interaction, the fermionic constraint takes a rather simple form: 

1 (jrg 



(2.3) 



While this is just a first order differential equation, the term coming from the interaction makes the 
analysis difficult. Still, in a classical level, this equation can be exactly solved [4], which helps to 
understand properties of the "chiral transformation" on the light front that is different from the ordi- 
nary one because the chiral transformation is imposed only on the dynamical ("good") component [5]. 
But, as was discussed in Ref. [4], in order to describe the symmetry breaking, we need to solve the 
fermionic constraint in a quantum level, which we will indeed do in the next section. It should be 
noticed that, in the quantum theory, the operator ordering must be specified, though it is a priori not 
known. Thus in Ref. [4], they defined the quantum theory by the above operator-ordering. Namely, 
V'+ comes to the left of color singlet operators such as or ip^ij^tp^. Since this operator ordering 

correctly reproduces physics of the chiral symmetry breaking, we adopt the same operator ordering in 
the present paper. 

With the vector interaction included, the corresponding fermionic constraint becomes terribly 
complicated, while its structure (as a differential equation) is the same as Eq. (2.3). We show its 
explicit form for later convenience: 



^ (i7±d±i + mo) 7'^V'+ - ^ 



+ Gy V-(: V'+7"^V'+ + 75V'-(: V'+7^75V'+ 



7+ ('(/'+?/'_ + - i757'^V'+ (■0+^5-0- + '0-n5'0+) 



+ 



717^-0+ : {^+7iV'- +V'-7iV'+} : "71757^^+ : {^+7175^- + ^-7175 V'+ 



Here we have followed the same operator ordering as in Eq. (2.3). This equation only makes sense in 
the quantum level because we have already taken the normal order. 

2.3 Quantization 

Since the bad component of the spinor is not an independent degree of freedom, the quantization 
condition is imposed only on the good component. From Dirac's procedure, one finds the following 
anti-commuting relations: 



{tp+a{x), i^ipiy)}^^^^^ = ^(A+)a/3 S{X - y) , 

where^ 5{x — y) = 5{x~ — y") 5^'^^ (xj, — y±). The mode expansion of the field is given by 
i^^{x) = I -^^-^Yl [b{p,s)u+ip,s)e-'P-^ + d\p,s)v+{p,s)e^P-^ 



(2.5) 
(2.6) 

(2.7) 



^Notice that x = {x ,xY) for space coordinates, but p = {p'^,p±) for momenta, and p ■ x = p a;"*" — px with 
px = —p'^x~ + kj_x''j_. See Appendix A. 
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where we have introduced notation for the integration 

/poo poo 
Jo J-oo 

and u^{p,s) and v+{p,s) are the good components of the solutions (with spin ,s) to the free Dirac 
equation, i.e., ?t+ = A^u and v^- = A^-V (for more details, see Appendix B). Notice that the lon- 
gitudinal momentum takes only positive values because of the special property of the dispersion 
relation on the light front p" = {p'j_ + m^)/2p+ > 0. The operators b{p, s) for quarks and d{p, s) for 
anti-quarks satisfy the anti-commutation relations: 

{bip,s), b\p',s')} = {d{p,s), d\p',s')}=6ss'Sip-p'), (2.8) 

and zeros for other combinations. The normal order in Eq. (2.4) is defined with respect to these 
operators. 

3 Fermionic constraints and LC Hamiltonians 

Here we solve the fermionic constraint in the quantum theory, and compute the LC Hamiltonian which 
is one of the necessary ingredients for deriving the bound-state equations. While the whole procedure, 
with the vector interaction included, is really tedious, it is straightforward and essentially the same as 
in the case without the vector interaction which was recently achieved in Ref. [4]. Thus, in this section, 
we first discuss the case without the vector interaction in order to demonstrate the methods we use. 
The presentation is basically along the previous analysis of Ref. [4] , but we here adopt the spinor basis 
for the mode expansion (2.7) (cf. instead of a simple Fourier expansion, Eq. (3.2) in Ref. [4]) and do 
not use any specific representations of the 7 matrices, both of which are much more convenient for 
the extension to the case with vector interaction. In order to avoid too much technicalities, but still 
to keep the paper self-contained as much as possible, we put some of the details aside to Appendix 
D. The involved case with vector interaction follows after this, with emphasis only on the differences 
coming from the inclusion of the vector interactions. Notice also that the analysis without the vector 
interaction can be applied to the scalar and pseudo-scalar mesons since we simply ignore the effects 
of vector interaction in these channels. 

3.1 The case without vector interaction as a warm-up 
3.1.1 Solving the fermionic constraint in bilocal form 

It was found in Ref. [4] that one can systematically solve the fermionic constraint by using the 1/N 
expansion after rewriting it with respect to color singlet bilocal operators. The 1/N expansion of the 
bilocal operator equations is made possible by use of the boson expansion method, which has been 
known as a standard technique in many body physics, but was only recently applied to the light-front 
field theories in Ref. [12]. The basic quantity is defined by a color singlet bilocal operator of the good 
component at different space points and at the same LC time: 

M„^(x,y) = V2VU(*)V'+/3(y), (3.1) 

where x = (x~,a;^) and a, (3 are the Dirac indices, and the color indices are suppressed. We further 
introduce bilocal operators of scalar (S) and pseudo-scalar (P) as follows: 

SK{x,y) = *(a;)A_*(y) + h.c, (3.2) 
P^{x,y) = *(a;)A_(z75)*(y) + h.c, (3.3) 
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where h.c. is the hermitian conjugate and the subscript "R" imphes that it is the "real" part (later 
we will define the "imaginary" part when we discuss the case with vector interaction). Notice that 
these bilocal operators reduce to the scalar and pseudo-scalar operators at the same space-time point: 
Sy{{x,x) = *(£c)^'(a;) and PR,(a;,a;) = ^ {x)ij5'^ {x) . Notice also that these quantities depend upon 
the bad component of the spinor ■0_ and thus are not known until we solve the fermionic constraint. 
In other words, once we rewrite the fermionic constraint (2.3) with respect to the bilocal operators, 
they can be interpreted as the constraint equations for these unknown bilocal operators. In particular, 
the equations for Sr and Pr form a closed set of equations. Indeed, one obtains in momentum space 

g+SR(p,g) = ^[YQi + mo)^^(Mai3{p,q)-Map{q,p)) 

7^ r^'^^^' 9 - - (Sr + Pr) Jk, I) 

- (Sr - i75 Pr)„^(^'0 Ma/3(Q -k- Z,p)| , (3.4) 
9^PR(p,g) = -^{(n5)(7'd + "io)}^^Ma/3(p,q) - \{{l\\ + mQ)(i-if,)^^^Uo,p{q,p) 

+ (i75SR + PR)^^(/e,Z)M„^(g-fc-Z,p)| , (3.5) 

where we have used short-hand notations like (Sr(p, q))^^ = (5q,/3 Sr(p, q), (Sr + ^75 Pr)q,^(/c, i) = 
(SR(fc, /) +^75 PR(fc, l))ai3, etc, and the momentum integration is from —00 to 00 (even for longitudinal 
momenta), which follows from our definition of the Fourier transformation 

Wc have written equations only for Sr and Pr (instead of a generic matrix with spinor indices con- 
structed by good and bad components like V'+aV'-/?); but these two equations are enough for the 
purpose of computing the Hamiltonian (see below). 

We will solve Eqs. (3.4) and (3.5) by using the 1/N expansion. We expand the bilocal operators 
as follows: 

MaMP>Q) = ^E(^j m(5(p,q), (3.6) 

and similarly for S{p,q) and P{p,q), and the expansion coefficients of them are written as s^'^\p,q) 
and p(")(p, q), respectively (for notational simplicity, we omit the subscript R in this subsection). 
Each contribution m^'^\p,q) is given by the boson expansion method of the Holstein-PrimakofF type 
[12, 4] as a function of a bilocal bosonic operator B{pi,si : P2,S2) satisfying 

[b{Pi,si : P2, S2), B'i{p[,s[ : P2, 4)] = ^si^ <5(Pi - p'l) Ss^s'^ KP2 - P'2) > (3-7) 

and so on. Note that there is no A''-dependence in m(")(p, q) since B{pi,si : P2, S2) is of 0{N^). Note 
also that this bilocal operator is in the spinor basis (i.e., with quantum numbers (spin) si, S2), and is 
not the same as the one in the Dirac basis (i.e., with the Dirac indices like a, (3) introduced in Ref. [4]. 
The spinor basis is more convenient for treating the vector mesons than the Dirac basis. The boson 
expansion method is constructed in such a way that the commutator which is satisfied by the bilocal 
operator Ma/3(p, g) is correctly fulfilled by using the bosonic operator B{p, s,p' , s'). It is useful to 
recognize that the commutator between two : M :'s becomes bosonic in the large N limit. Thus, in this 
limit, the bosonic operator B (or B^) becomes identical to the bilocal operator : M : (depending on 



8 



the sign of the momenta) and thus can be understood as annihilation/ creation operators of a quark- 
antiquark pair. The definition of this expansion and exphcit form of the first few terms are given in 
Appendix C. 

Substituting the lowest order expression (C.ll) of the boson expansion, one can compactly rewrite 
the lowest order 0{N) of the fermionic constraints (3.4) and (3.5) as 

where we have rescaled the coupling constant Gs = 5s /-^ (Recall that the large N limit should be 
taken with GsN being fixed). The solutions are easily obtained: 

s(°)(p,qr) = -M'-^5{p + q), (3.9) 
p(°)(p,qr) = 0, (3.10) 
where M is the dynamical mass of the fermion and is given self-consistently by 



M = mo - Gs (0| |0 
d^p d?q 



"^'-'^I sW(p,g) + 0(l/ViV). (3.11) 



As we will see below, this equation becomes the gap equation after we impose an appropriate cutoff. 
Substituting the solution (3.9) into the above definition of M, one encounters an integral which is 
divergent both in infra-red and ultra-violet regimes: 

M-mo_ f d^p e(p+) 

9s / TTT-Ts— ^ • I3.i2j 



M (27r)3 p 

We regularize this integral by introducing a cutoff A. We have to admit that even physical quantities 
may have explicit dependence upon the cutoff A because the NJL model is not a renormalizable theory. 
Besides, it is important to know that the cutoff cannot be taken arbitrary. The key for finding an 
appropriate cutoff is to keep as many symmetries of the system as possible. It turned out [4] that the 
parity^ invariant cutoff 1^=*= | < A yields a reasonable result (even in the chiral limit) 

where we have used the dispersion relation of a free fermion with the dynamical mass = 2p'^p~ — 
p\ = to derive the lower limit of the p'^ integral. Imposing a cutoff on both the LC energy and 
the longitudinal momentum seems to violate the boost invariance of the system, but this is not the 
case because is an internal momentum to be integrated out. We will, however, have to be more 
careful when we apply the similar cutoff to the two body sector. 

When the coupling constant is larger than the critical value ^g'^^'*^ = 27r^/A^, the gap equation 
(3.13) develops a nonzero solution even in the chiral limit. This is easily confirmed by the numerical 
solution to the gap equation, as is shown in Fig. 1. Therefore, if we take this nonzero solution^ and 
substitute it to the canonical Hamiltonian, we are able to describe the dynamics in the broken phase. 

®The parity transformation in the z direction corresponds to the exchange <-> p~ in momentum space. 
Strictly speaking, the symmetric solution for > g'g'^''*^ appears only in the chiral limit. If the current quaxk mass 
is taken non-zero and is decreased to zero, then the solution naturally approaches to the broken solution for > Pg"'*\ 
and to the symmetric one for gfg < gg'^'^''^ This is evident on Fig. 1. The symmetric solution for stronger coupling 
is unstable against the small explicit violation of the chiral symmetry, just like in the Ising model under the external 
magnetic field. 
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M/A 




Figure 1: Dynamical quark mass M/A as a function of the dimensionless coupling constant (7sA^/47r^ 
for different current quark masses, mo/A = 0.08, 0.06, 0.04, 0.02, (chiral limit), from top to bottom. 

As we will discuss shortly, solving the bilocal fermionic constraints (3.4) and (3.5) up to the next- 
to-leading order is important for deriving the LC Hamiltonian to the first nontrivial order. In fact, one 
can even solve the constraints at arbitrarily high order, in the inductive way. Namely, if one knows 
the solutions up to the n-th order, then the {n + l)-th order solution is given by some known functions 
and the solutions whose orders are lower than n + 1. All the details are discussed in Appendix D, 
together with the derivation of the LC Hamiltonian. 

3.1.2 The LC Hamiltonian 

Since the canonical LC Hamiltonian is bilinear with respect to the spinor (see Eq. (D.5)), one can 
rewrite it by the bilocal operators: 

= IJd'pJ d\ {iq'MiP, q) 6{p + q) + '^ j d'p Sr(p, -p) , (3.14) 
where we have introduced a new bilocal operator, 

iyiix,y) = ^(x) A^-f{^{y) -h.c. . (3.15) 

Note that this new bilocal operator can be related to the known operators, M, S, P (see Appendix D). 
Therefore, given the solutions and \ one can immediately obtain Vj ^"'\ the expansion coefficients 
of Vp Hence, one can compute the LC Hamiltonian order by order 

^-~|(^)"4"'. (3.10, 

where the subscript S was put in order to remind that this is for the case with Cs (2.1) (without the 
vector interaction). We substitute the "broken" solutions of Sj^^ and Pj^^ to obtain the Hamiltonian 
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which describes the broken phase of the chiral symmetry. It turns out that the lowest order h^^^ is 
divergent but is just a constant, and thus can be neglected for the present purpose. The next order 
vj {p, —p) and s^^-* (p, —p) become strictly zero, and we have h!^^ = 0. Therefore, the first nontrivial 
contribution is given by the next-to-leading order n = 2. After straightforward calculation, one arrives 
at the following expression (see Appendix D for the derivation): 



r f „2 , ^2 -I- M2 1 

J [dp] [dq] I ^^^-T- + 2q+ \ ^ ^' ' ^' 

f [dp] [dq] [dk] [dl] Sjp + q-k-l) ^ ^ . , , 

|?2(p, si)i;(q, S2)} [v{l, S2)u{k, 4)} + {?^(p, si) ^75 ^^(g, S2)} {^^(i, 4) ^75 ^^(fc, 4)} 



,(3.17) 



where we have introduced a constant 



As already mentioned, in the leading order of the 1/A^ expansion, the bosonic operator B^{p, si : 
q, S2) can be interpreted as a creation operator of a quark-antiquark pair (the pair behaves like a 
boson). Since this Hamiltonian contains only terms of B^B type, it describes the dynamics of the two 
body qq sector. Indeed, the first, diagonal, term corresponds to the free part for this quark-antiquark 
pair with constituent mass M, which is easily understood by the fact that {p'^ + Af^)/2j»^ + (g^ -|- 
M^)/2g+ corresponds to the sum of LC energies of a free quark and a free antiquark. The second term, 
which is proportional to the (scaled) coupling constant and thus clearly comes from the interaction, 
needs to be diagonalized. As we will see later, diagonalization leads to a bound-state equation for the 
quark and anti-quark pair. 

It is of special importance to recognize that the divergent factor Kg defined by Eq. (3.18) is 
independent of the external momentum P+. This is formally seen by changing the variable in the 
longitudinal integration (y = r'^/P'^) 

dr-^^=rdy'^'^ 



P+-r+ i-00 1 - y ■ 

Of course we need to regularize this integral, but if we use a cutoff which does not depend on the 

external momentum , then the result docs not depend on P"*", either. This is true as far as wc treat 
a boost invariant cutoff since the external momentum P~^ = p'^ + q~^ is eventually related to the total 
momentum of the two body qq state (see discussion in the next section). 



3.2 The case with vector interaction 

Inclusion of the vector interaction certainly complicates the fermionic constraint, but the structure of 
the constraint (as a differential equation) is the same and we can follow the same procedure as before. 
Thus, we do not repeat the whole procedure here, but rather present only the differences from the 
previous case. The main difference is that we have to treat another types of bilocal operators, which 
makes the number of equations double. 



3.2.1 Definition of the bilocal fields 

In the case without the vector interaction, the fermionic constraint (2.3) was transformed into two 
equations (3.4) and (3.5) which form a closed set for bilocal operators Sr and Pr. On the other hand. 
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when the vector interaction is turned on, the original fermionic constraint (2.4) cannot be expressed 
only with respect to Sr and Pr, but we need to introduce another types of operators, the vector 
and the axial vector A^^. Specifically, we define the following bilocal operators: 

Six,y) = ^{x)A_^{y), P{x,y) = ^{x)A_iij5)^{y), (3.19) 

y\x,y) = ^ix)A_{ji)^{y), y^ix,y) = h{xh^^{y), (3.20) 

A\x,y) = t(a;)A_(7i75)*(y), AHx,y) = h{x)j^j5^iy). (3-21) 
For each operator = S, P, V^, A'^, we further define its "real" and "imaginary" parts as follows: 

0{x,y) = ^{OKix,y) + iOi{x,y)}, (3.22) 

where 

OK{x,y) = Oix,y) + Oix,y)\ 
iOi{x,y) = 0{x,y) -0{x,y)\ 

so that they satisfy OR,^i(a;, y)^ = ORj(a;, y ). Of course these definitions are consistent with (3.2), (3.3) 
and (3.15), which were already introduced when we solved the fermionic constraint without vector 
interaction. Note that all the "real" parts at the same space-time points (i.e., OR(a;,a;)) reduce to the 
operators which have appropriate statistics. Namely, 

SR(a;, x) = ^{x)^{x), Pn{x, x) = ^{x) ij^ ^{x), (3.23) 
y^{x, x) = *(a;) 7^ *(a;), A^(a;, a;) = ^(a;) 7''75 ^{x). (3.24) 

Actually, the operator 0{x,y) has been constructed so that Or(£C,£c) becomes the standard bilinear 
operator with appropriate statistics. These "interpolating" operators become useful later when we 
determine the kinematical structure of the lowest Fock state of a meson. 



3.2.2 The bilocal fermionic constraints 



Having defined the necessary bilocal operators, we are able to write down a closed set of the bilocal 
fermionic constraints, which are obviously much involved than the previous case. A closed set is 
formed by S, P and the transverse components of vectors V* and A*. Indeed, the explicit form of the 
bilocal fermionic constraints are given by 



S(p,q) 



A\p,q) 



'^^q'^ + mo 



af3 



Mal3iP,q) 



+ 



- Ma/3(P,Q - k-l) 



1 



a/3 



{k,l) 



, (3.25) 



-7*75 
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+ 



. (3.26) 



The plus components of the vectors V+, A+ do not contain the bad spinor and are expressed by M„^, 
while the minus components V~, A~ are expressed in terms of S, P, V* and A*. We can solve the 
above equations following the same procedure as before, though all the calculations arc tremendously 
cumbersome. It should be noticed that the leading order constraint for S again leads to the gap 
equation for the dynamical quark mass M, which is equivalent to Eq. (3.13) because we took the 
normal order in the definition of the vector interaction (2.2). 

3.2.3 The LC Hamiltonian 

Substituting the solutions of the bilocal fermionic constraints into the canonical LC Hamiltonian, 
one finds again that the first nontrivial Hamiltonian starts from n = 2 in the 1/A^ expansion (see, 
Eq. (3.16)). As a result, the Hamiltonian gets extra terms due to the vector interaction in addition 
to /ig given by (3.17). After a lengthy calculation, one finally obtains (unimportant c- number is 
ignored) 

(2) ^ {dv\\dq\\dk\\dl\ 5ip + q-k-l) 



(27r)3 4^/p+q+k+l+ 
H ('^v + '^SV + + ^sy^)^Hp^ Si ■ 9> S2)B{k, s[ : I, 4) ■ (3-27) 



Sl,S2 .s[,s!2 



The "kernel" cj's in the interaction terms are functions of momenta p,q,k,l and spins si, S2, s'l, S2 

and given by 

uj^j = gylny Y. + Y \[{u{P,sih''v{q,S2)}[v{l,s'2)j;j,u{k,s[)} 



+ 



{^Z(p,si)7'^75i;(g,S2)}{i)(Z,4)7/.75w(fe,s'i)}] , (3.28) 



^SV = Ws9yKs 



dbdc e{b+)M 



_ 2 
= Qy 



(27r)3 6+C+ 

5{p + q + b + c) ^u{p, si)7+75t;(q, S2)}{u(Z, s'2)i^^u{k, si)| 

+^(fc + ^-6-c)|^Z(p,Sl)^75^;(Q,S2)}{^'(^,S2)7^75^i(fe,s^)} ^ (3-29) 
[ dbdc e(6+)(62 +m2) , ^ 

{u{p, si)7+f (q, S2)}{v(Z, s'2)-i^u{k, s^)} + {u(p, si)7^75^^(g, S2)}{v{l, S2)7^75^i(fe, s'l)}] , 
dbdcdb'dc' e{b+)M e{b'+)M 

(3.31) 



x5{p + q + b + c) 5{k + I ~ y - c!)^u{p, si)7+75i;(q, S2)^\v{l, S2)7"^75«(fe, s'l)} , 
where we have introduced another type of divergent constant 

d^r e(r+) 



{l+»v/ 



-1 



(27r)3 P+ - r+ 



(3.32) 



Notice that Ky is independent of the external momentum P+ due to exactly the same reason as for 
Kg given by Eq. (3.18). 
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A few comments are in order about this LC Hamiltonian. First of all, we have introduced the 
rescaled vector coupling constant = GyN which naturally yields the expansion of the LC 
Hamiltonian. Otherwise all the interaction terms become larger than h'g^ in the 1/N counting, which 
is not interesting. Next, it is clear that all the interaction terms are of the B^B type. This means 
that the Hamiltonian does not mix the Fock components with different number of particles. Namely, 
we can solve the Hamiltonian within a small Fock space with one boson (or, equivalently, a qq pair). 

Lastly, it would be necessary to explain why we obtained terms like wsv and Wgv^ which are of 
the higher order with respect to Gs and Gy, which seems less trivial compared to the case without the 
vector interaction. The mathematical reasons are the following: First, the bilocal fermionic constraints 
(3.25) and (3.26) have contributions proportional to Gs or Gy. Thus, it is natural that the solutions 
to the fermionic constraints, in each order of the 1/N expansion, can depend upon Gs or Gy- On 
the other hand, the LC Hamiltonian contains Vf, and this vector is again represented by other bilocal 
operators. Therefore, if we substitute the solutions (which depends on Gs and Gy) to Vj in the LC 
Hamiltonian, there appear higher-order terms with respect to Gs or Gy. 



4 Bound-state equations of light mesons 

4.1 Interpolating fields and the lowest Fock state representation of mesons 

In the leading order of the 1/N expansion, mesonic states arc expressed as "constituent" states with 
a dynamical quark-antiquark pair, as was discussed in Rcf. [4]. For the purpose of constructing a 
meson with appropriate statistics, the first thing to do is to determine its kinematical structure. This 
can be done in such a way that the mesonic state has the same structure as that of its interpolating 
field. Namely, a part of the LC wavefunction which rather trivially depends on the statistics can 
be determined even withoTit solving the dynamical bound-state equation. Then, after extracting this 
known part, we can define the spin-independent part of the LC wavefunctions. 

As we already discussed, the interpolating fields are given by the "real" parts of the bilocal op- 
erators at the same space-time points (see Eqs. (3.23), (3.24)). We can find their explicit forms in 
terms of spinors and the bilocal boson operators, by using the solutions to the bilocal fermionic con- 
straints. Since (in this subsection) we are interested in the kinematical structures which is independent 
of the dynamics, we are allowed to use the solutions to the fermionic constraint without the vector 
interaction^. To the leading order, one finds for the scalar and pseudo-scalar fields 



#(£c)z75*(a;) 



SR(a;,a;) 
Pr{x,x) 



[dP] 



iPx 



d^k\ 



2^k+{P-k)+ 

J2 {uik,si) (^.l^^^viP-k,S2)}B^{k,si:P-k,S2) + h.c. 



Sl,S2 



(4.1) 

+ 0(iV°), 



and for the vector 



^{x)-f^'^{x) = y'^{x,x) 



^ [u{k, si)-ft'v{P - k, S2)}B\k, si -.P-k, S2) + h.c. 



Sl,S2 



(4.2) 



+ 0{N°). 



In fact, we can even use the free theory (a free fermion with dynamical mass M) to determine the kinematical 
structure. However, we do not perform this because we already have the solutions to the case with Cs =^ 0. 
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The four interpolating fields of the vector are obtained separately for the LC components ^ = ± and 
the transverse components n = 1,2, but it turns out that they can be combined to give the reasonable 
representation as above. The mismatch in the expressions above between (S, P) and V (namely, a 
factor Kg in Eq. (4.1)) is due to our use of the solutions with Cs and Cy = 0, and would be absent 
if one used the free field solutions. Needless to say, this difference is irrelevant for the present purpose 
of determining the kinematical structure. 

Since the interpolating fields are all defined at a single spatial point, they can be interpreted as a 
kind of point-like states without spatial structure. However, what we are interested in is the spatial 
structure of mesons, and it is indeed the (spin-independent part of the) LC wavefunction which is 
responsible for this spatial structure. (In other words, for the interpolating fields without spatial 
structure, the LC wavefunctions are simply constant in momentum space.) First of all, the mesonic 
states having total momentum P can be represented as follows {x = k'^/P'^): 

I meson ; P ) = j dx I d'^k± \ ^meson 

(x, k±, si, S2)B^{k, si:P-k, S2) |0) , (4.3) 

where $meson(a^! k±, si, S2) is the LC wavefunction of the qq Fock state, and the vacuum is defined by 
the bilocal boson operator 

B{p,Si:q,S2)\0)=0. (4.4) 

Notice that this vacuum is equivalent to the one defined by the quark antiquark operators to the 
leading order of the 1/A^ expansion. From the interpolating fields shown above, we can determine the 
kinematical structure (or, spin-dependent parts) of the LC wavefunctions: 

^T,{x,k±,si,S2) = -—^^==(j)n{x,k±)\u{k,si)i-f5v{P-k,S2)], (4.5) 
2^/x[l — X) ^ ^ 

^aix,k±,si,S2) = ^ 1 } .(^a{x,kx_)[u{k,sx)v{P -k,S2)\, (4.6) 
2-\/x(l — x) ^ ' 

^^p{x,k^,si,S2) = ^ . } 6p{x,ki_)e^{\,P){u{k,si)-f''v{P-k,S2)], (4.7) 

2yjx(\ — X) ^ J 

where we have introduced the spin-independent LC wavefunction^ 4'{x,k^) for each meson. These 
expressions are valid when P± = . The corresponding formulae for P± ^ are shown in Appendix E. 
For the vector meson, we have already introduced 3 physical states \p,\) with helicity A = ±1, (see 
Eq. (4.7)) which should be distinguished from 4 states \p^) naively constructed from the interpolating 
field without the polarization vector e^(A, P). Below, we will explain the necessity of using |p. A) 
instead of \p^)- 

One can directly verify the orthogonality among these mesonic states. For example, let us check 
the orthogonality between the scalar \a) and the vector \ff') when the transverse momentum of the 
meson P± is zero: 

(c7;P+,Px = 0| ff;P'^,P'^ = Q) 

= 5{P-P')\p^=^P+ j^^d^k^<Pl{x,k^)<Pp{x,k^) 

{u{k,si)rv{P-k,S2)}{v{P-k,S2)u{k,si)}\^ (4.8) 



X 

Sl,S2 



where the quantity on the last line can be evaluated as 

J2 si) 7^ v{P - k, S2)}{v{P - k, S2) u{k, Si)} I 



S1,S2 



P±=0 



The spin-independent LC wavefunction in the present paper is different from that of the previous papers [4, 5] by a 
factor 2x{l — x). 
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4MP+(l-2x), ^ = + (4.9) 

AM (l-2a:) 
2PTx-(l-x) 



Assuming that the spin-independent wavefunctions arc functions of x(l — x) and fe^ (which can be 
checked a posteriori, but is natural due to the charge symmetry), one concludes that the bracket 
{a\ p'^) above is zero for all /x. The orthogonality for other combinations of the states are verified too 
in a similar way. 

It should be noticed that the four vector states \p'^) are not independent of each other (thus 
not necessarily orthogonal to each other), since a real massive vector meson should have only three 
physical modes, i.e., two transverse and one longitudinal modes. This can be seen by the fact that 
the spinor part of the state (4.7) satisfies the following identity: 

{Ek + E^p^k)) Hk, si)7+ v{P - k, S2)] + P+ [u(fe, si)j- v{P - k, ss)] 



-P 



± 



u{k,si)Y v{P - k,S2) 



= 



where E^ and Ep^k are the LC energies of the quark and the antiquark 

The projection onto three physical modes can be easily done by the circular polarization vector 
e^{\,P). As we already introduced in Eq. (4.7), three independent states with helicity A = ±1,0 
are given by 

|p,A;P) = e^(A,P) |p^P). (4.11) 

These three states are orthogonal: ( />, A | p, A' ) = for A ^ A'. In the frame where the transverse 
momentum of the vector meson is zero, P = {P^,P±) = (P~^,0±), the polarization vectors are given 
by(e'^ = (e+,e-,ei,e2)) [13] 

e'^(A = ±l,P) = (0,0,T-^,--J=), (4.12) 
e''{X = 0,P) = (£^,-^,0,0), (4.13) 

satisfying e^(A, P)e^(A', P)* = —Sw, P^e'*(A, P) = 0. Therefore, the explicit relations between \p>^) 
and the physical states \p, A) are 

|ftA = ±l) = ±-^=|p') + ^l(>'>, (4.14) 

k,A = 0) = Z^|,+ ) + ^|,-). (4.15) 

Notice that these relations are valid only in the frame P = (P+,0_l), and does not hold any longer if 
one moves to other reference frame. 

In summary, the orthogonality holds among the five mesonic states: |o") , |7r) and |p, A) with 
A = ±1,0. Therefore, together with the normalization of each state, we obtain 

{a,P \a',P') = 16Tr'^P+6a,a'SiP-P'), (4.16) 

where a {a') distinguishes the species of meson a = cr, tt, (p, A). This also leads to the normalization 
condition for the LC wavefunctions: 

1 rl roo 

— dx d^^i ^ |$a(x,A;x,si,S2)P = l. (4.17) 

iOTT Jo J-oo 



S1,S2 



16 



4.2 Bound-state equations 

Before going into details, let us briefly discuss a generic qq state in order to clarify the procedure we 
perform. In general, the LF energy of the two body state may be schematically written as 



2k+ ^ 2(P+ - k+) 



P,, = -^^TT^ + ' tr.^ ' + P), (4.18) 



where the first two terms are the "kinetic" energies of the quark and the antiquark, and V is the 
potential which allows for a bound state. This form of the energy leads to the following bound-state 
equation: 

f 2 kl + M^' 



^meson 



x{l — x) 



(p{x,k±) = J^ dyJ(fp±V{x,k±;y,p±)(l){y,p±), (4.19) 



where we have chosen P = (P+,Pj^) = (P"^, 0_l) for simplicity, and redefined V with some factors 
included. 

For the scalar and pseudo-scalar mesons, the potential V has a simple structure, and we can 
evaluate it without any further approximation. On the other hand, for the vector mesons, V has 
nontrivial dependencies upon g^. For the present purpose of seeing the leading effects of the vector 
interaction (2.2) on the vector channel, it is actually enough to derive the potential V up to the 
leading order of the vector interaction g^. We will see that in this leading order the potential term 
V{x,k±;y,p±) is separable with respect to the internal iy,p±) and external {x,k±) variables, and 
actually depends only on y and p±. 

In the following, we are going to derive a bound-state equation for each meson in the frame 
P = (P"*", 0_l), which is simply deduced from the light-front eigen- value equation: 

2 

P- I meson ;P+,Px = 0) = | meson ; P+, P± = O) . (4.20) 

To the first non-trivial order in the 1/A^ expansion, this equation can be restricted to the qq subspace 
of the Fock space, leading to the bound-state equation for mesons. 

4.2.1 Scalar and pseudo-scalar mesons 

It is straightforward to derive the bound-state equations for scalar and pseudo-scalar mesons. In 
Eq. (4.20), we use the LC Hamiltonian P~ = h^^ defined in Eq. (3.17) (recall our strategy for scalar 
and pseudo-scalar mesons, see Sect. 2.1) and the scalar and pseudo-scalar states given by Eqs. (4.6), 
(4.5). Then, one obtains 



mi - 



A;2 + m2 



^(->fc-) = J ^V "';i\"_'f^% .(y,P.), (4.21) 



xil-x) j^"' ' (27r)3yo vHl-y) 



An important property to be noticed is the fact that, in both cases, the potential terms (the right-hand 
sides) are independent of the external variables x, k±. This property is due to the simple structure 
of the contact interaction £g. The same thing will be seen in the case of the vector meson. These 
results are equivalent to the previous ones obtained in Ref. [4]. Formal solutions to these equations 
are easily obtained (due to the simple structure of the potential terms), but we postpone to discuss 
them until the next section, where we present the LC wavefunctions and masses of the mesons after 
careful evaluation of divergences. 
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4.2.2 Vector mesons 



For a vector meson, the first nontrivial contribution (in the 1 /N expansion) of the LC Hamiltonian is 
given by P" = h^"^^ in Eq. (3.27), and the transverse and longitudinal states are given by Eqs. (4.14) 
and (4.15), respectively. Although the mass and the spin-independent LC wavefunction must be, 
by definition, the same for both the transverse and longitudinal polarizations, we distinguish the 
bound-state equations for these two modes since they appear differently until we perform careful 
regularization of the divergences. Therefore, in the following, we use notations like Vx/L) 0t/l(^) ^±)j 
and TTT-T /L for corresponding transverse and longitudinal quantities (though eventually they are proven 
to be equivalent to each other). 

For a transversely polarized vector meson, a lengthy calculation yields the following potential Vr 
(in the sense of Eq. (4.19)) 



Ft = - 



(27r)' 



1 + 



(27r)^ 



d\± 



P+ 



1+ 



pI + M^- 2y{l - y)p'\ 

y2(l _ yY 



(4.23) 



This potential comes from tOy in the LC Hamiltonian (3.27). The other terms do not contribute to the 
bound-state equation because of the orthogonality among the bispinors. Notice that the potential is 
already independent of the external variables x,k±. Taking the leading contribution of gy, one arrives 
at an equation for the LC wavefunction (p^^ix, k±) = (f>\=±i{x, k±): 



m.'Y — 



kl + M^' 
x{l — x) 



(f)rj^{x,k±) = 



_9w_ 

(2vr)2 



,2 pi + M 
p± — 



2y(l - y)p] 



y2(l - y)2 



My,p±)- (4.24) 



Since we pick up the leading contribution in Eq. (4.23), the divergent factor in the square bracket is 
now replaced by 1. This factor is the explicit expression of /Cy which appeared in the Hamiltonian 
(cf: see Eq. (3.32)), and if we keep this factor, the structure of the bound-state equation looks similar 
to those of the scalar and pseudo-scalar mesons (their potential terms are proportional to g^Ko^)- 
However, as we will see below, the corresponding equation for the longitudinal mode does not have 
the same simple structure, and we keep the leading order effect (namely, Ky = 1) for simplicity. 

The longitudinally polarized state is much more involved. One of the main complexity comes from 
the fact that V"*" and V~ are treated (and in fact do behave) differently in the light-front quantization, 
while the physical longitudinal mode is a linear combination between these two. This situation requires 
a longer and attentive, but still straightforward, calculation which eventually leads to the following 
potential Vl (for the details of the derivation, see Appendix F): 



gy 
(27r)3 

2 + 



gy 



4{kl + M^) r gy /-i /• 2 



ml x{l — x) — {kj_ + . 




2 fe^ +M2' 




-1 



It should be noticed that the interaction term iOy^ in the LC Hamiltonian (3.27) as well as Uy 
contributes to the longitudinal potential while it does not contribute to the transverse potential. As 
we mentioned above, the gy dependence is not the same as the transverse case, and is much more 
non-trivial. However, we obtain a simpler result by taking the "leading term" with respect to gy in 
this potential (more precisely, in the eigen- value equation). Here we give a quick, but less systematic, 
derivation of the final result. (A more detailed derivation is explained in Appendix F.) This can be 
done as follows: First of all, we ignore the gy dependent term in the second line of Eq. (4.25) that 
gives the higher order in ^y and thus can be ignored anyway. We immediately find 



ki 



x{l — x) 



' k^+M^' 



x{l — x) 



4 (fc2 + M2 



m5 



x{l-x)-{k'i + M'^) 
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Then we integrate the resulting bound-state equation over x and k±, obtaining the fohowing: 



J dx(fk± l^ml - ^"^j^"^-^ I (Phix, k±) = "4^^^ (^J dxd^k±^ J dy d^p± ^^'^ '/'l ^ P± ) • 

Here we have kept the gy dependent term in the curly bracket in the first line of Eq. (4.25). We can 

modify the bound-state equation by using this relation in the right-hand side of it. Taking the leading 
term with respect to ^y, we finally obtain the following simpler equation for the longitudinal mode 
^Lix,k±) = (l)\=o{x,k±): 

It is evident that the right-hand side is again independent of the variables x,k±. The derivation of 
Eq. (4.26) looks tricky, but we can obtain the same result by carefully taking the leading order effect 
with respect to gy in the eigen- value equation (see Appendix F). 

Now we obtained the bound-state equations for the transverse and longitudinal modes (i.e., 
Eqs. (4.24) and (4.26) ). However, at first glance, these two look different from each other and 
thus seem to give different masses for the transverse and longitudinal vector mesons. This is of course 
physically unacceptable, and as we will verify soon, these equations are in fact identical and give the 
same mass rrirp = ruj^. This equivalence will be achieved after one specifies cutoff scheme. It is not 
hard to identify the origin of this (fake) discrepancy with the lack of Lorcntz covariancc in the light- 
front formalism. Namely, the main source is the fact that the LC components and the transverse 
components V* seem to behave differently in the light-front quantization. 



5 LC wavefunctions and masses of light mesons 

5.1 General arguments: Eigenvalue equation and the cutoff scheme 

In the previous section, we have found that all the bound-state equations have the potential terms 
which do not depend upon the external variables x, k_i_}^ This property leads to a significant conse- 
quence that the spin- independent LC wavefunction (f){x, k±) should have the same form as a function 
of x,k±. Let us see this in general arguments. Consider Eq. (4.19) with such a "separable" potential 
(the subscript i stands for cr, tt, T, L): 



x{l — x) 



(pi{x,k±) = dy J d'^p±Vi{y,p±)(j)i{y,p±) 



Since the right-hand side is just a constant, one immediately finds that the spin-independent LC 
wavefunctions should have the common functional dependence upon x and k± as we announced in the 
introduction: 

(J. 

'^i x{l-x) 

The constant factor Ci simply represents the right-hand side of the bound-state equation, 

Ci = J^dyJ d^p^ Viiy,p^_) (l)i{y,pi_) , (5.2) 

and is determined from the normalization condition of the wavefunction (4.17). Note that our LC 
wavefunction (5.1) is the simplest wavefunction from the general consideration of the relativistic 

^"This will be due to the contact interaction in the NJL model, but actually, within the present formalism, it appears 
in a highly nontrivial way, especially for the vector mesons. 
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bound-state equation [14, 15]; All the LC wavefunctions must have the similar denominator as in 
Eq. (5.1), while the numerator could be a function of x or k± in general. The spin-independent LC 
wavefunction always has a peak at a; = 1/2, but its shape (width of the peak) will change depending 
on the value of mass mj. Namely, the peak becomes sharp as nii becomes large. 

Substituting Eq. (5.1) into the above definition of Cj, one obtains an equation for (the eigen- 
value equation): 

2/(1-!/) 



2 rA 
mf ^ 



Though all the potential terms Vi {i = a, tt, T, L) are always negative, this equation makes sense 
because the denominator can be either positive or negative. 

Since the integral (5.3) is divergent in general, we need to regularize it. We evaluate the integral by 
introducing the "extended parity invariant cutoff' [4, 10] which is actually equivalent to the Lepage- 
Brodsky cutoff [16]: 

Indeed, this is a natural extension of the parity invariant cutoff in the two body sector, K~^K < 
where and K~ are the sum of (on-shell) quark and antiquark longitudinal momenta and energies 
[K+ =p+ + (P+ -p+) = P+, K- = {pl + M^)/2p+ + {p\ + M^)/2{P+ -p+) ]. This cutoff obviously 
preserves transverse rotation and parity symmetry separately, but in fact it does work better. First, 
it also respects the usual three dimensional space rotation [17]. Thus one can relate the above cutoff 
A to the 3-momentum cutoff I]i=i, 2,3(^0^ < ^sm through 2(A|i^ + M^) = A^. Next, the cutoff 
K^K~ < A^ is invariant under the boost transformation —>■ e^^K^, which is necessary for the 
relativistic formulation Therefore, we are going to evaluate the following regularized eigenvalue 
equation: 



Jy- Jo 



2/(1-2/) 



where 



R- 1 



Below, we will see that this cutoff derives reasonable results for the masses of mesons. 



5.2 Scaleir and pseudo-scaleir mesons 

We evaluate the eigenvalue equation (5.5) for the scalar and pseudo-scalar potentials: 

y ^ gs^S + (1 - 2j/)'M2 ^ ggKg pI + M^ 

(27r)3 y2(i_y)2 ' - (2vr)3y2(i_y)2- 

When we treat Eq. (5.5), we encounter two types of divergence in the longitudinal integrals. One is of 
the type dy/y which appears as the overall y integral in Eq. (5.5), and the other is dy/y from 
the factor Kg in the potentials. The first type of the divergence is regulated by using the extended 
parity invariant cutoff as indicated in Eq. (5.5). On the other hand, we already know how to regulate 
the second type of the divergence, since it appeared in the gap equation (3.12). However, in fact, we 

^^If one imposed the same condition \K^ \ < A as the one for the gap equation, it would violate the boost invariance, 
though yielding the same Lepage-Brodsky cutoff. However, as we mentioned before, there is no problem in putting 
I < A in the gap equation because the momentum is not the external momentum but the internal one to be integrated 
out. 
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do not have to know the exphcit expression of the regulated integral (regulated by the parity invariant 
cutoff) for the present purpose. We rather utilize the gap equation (3.12) to replace the divergent 
integral by a finite quantity. Indeed, the divergent factor Kg is modified so that it contains only the 
first type of the divergence: 
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(5.7) 



where we have used the gap equation (3.12) in the last line. For the remaining longitudinal integral, we 
can use the same cutoff as indicated in Eq. (5.5). For example, the equation for the pion is evaluated 
as follows: 
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In a similar way, the equation for nia is simplified to 
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(5.8) 



(5.9) 



These integral equations (5.8) and (5.9) coincide with the results obtained in the previous paper 
(cf: Eqs. (4.23) and (4.24) in Ref. [4]). The integral is analytically doable with the cutoff indicated 
in Eq. (5.5). Hence, we arrive at the following non-linear equations for the masses of scalar and 
pseudo-scalar mesons: 
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(5.10) 
(5.11) 



where we have introduced (square of) the ratios of the meson mass to the threshold mass 2M: 

2 / ^ \ 2 



\2Mj 



\2m) 



(5.12) 



Let us first look at the chiral limit mp = 0. Then, one finds trivial solutions = 0, ro- = 1. Namely, 
the pseudo-scalar meson becomes massless (the Nambu-Goldstone boson), while the scalar meson 
appears as a loosely bound state = 2M. 

When mo ^ 0, the equation for the scalar meson does not have a (real) solution for To- < 1 because 
the quantity in the curly brackets cannot be negative for < I- Therefore, the scalar meson does 
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Figure 2: Pseudo-scalar meson mass squared {m^/ K)^ as a function of mo/A. The solid line is the 
numerical solution to Eq. (5.10) for M/A = 0.4, and the dashed line is its approximate solution, 
Eq. (5.13). Two dotted lines are the threshold line (2M/A)2 = 0.64 (horizontal) and the upper limit 
of the bound state (mo/A)upper = 0.137- • • (vertical) as given in Eq. (5.14). 



not appear as a bound state^^. This observation is consistent with the literature [18, 9]. On the other 
hand, the pseudo-scalar meson is still a bound state for small nonzero mo and its mass is estimated as 



N rriQ 



1 



N 



1 1 1 + 



47r2 



(5.13) 



In order to compute the decay constant /t^, one has to determine the overall factor C^r of the LC 
wavefunction (5.1) from Eq. (4.17). Then, one finds /tt = 2MZ7r , which, combined with Eq. (5.13), 
implies the Gell-Mann, Oakes, Renner relation m^/^ = — 4mo(^'^'), as was discussed in Ref. [4]. In 
Fig. 2, we show the numerical solution to the eigenvalue equation (5.10), together with the approximate 
solution (5.13). Shown is the square of the mass, and thus the linear dependence upon mo is clear for 
small mo. As mo is increased, the pseudo-scalar meson's mass increases, and reaches at the threshold 
mTT = 2M. The value of bare mass mo which gives the threshold mass is easily calculated from 
Eq. (5.10) where we use = 1 and delete dependence by using the gap equation. Explicitly, the 
result is 

'mo\ M \^ 2- (M/A)2(l + ln■ 

X 



upper 




2(M/A)2lni±| 




(5.14) 



upper 



0.137- 



which is in 



For M/A = 0.4, this upper bound of the bare mass becomes (mo/ A) 
agreement with the numerical result. Therefore, the bound state in the pseudo-scalar channel exists 
even for nonzero bare mass mo if it is smaller than (mo/A)upper given by Eq. (5.14). 

Once we obtain the mass of mesons, we can discuss the shape of the spin-independent LC wavefunc- 
tions. It is the meson's mass that controls the shape of wavefunction. Consider the chiral limit, only 



^^The presentation in the previous paper [4] was misleading about this point. There, it was argued as if the scalar 
meson appeared as a bound state even for mo 7^ case. It does not make sense to write — 4M^ + ©(mo) (i.e., 
Eq. (4.29) in Ref. [4]) for mo 7^ within the present approximation. 



22 



[X, kh 



1 




0.8 



0.6 



0.4 



0.2 



0.2 



0.4 



0.6 



0.8 



1 



X 



Figure 3: The spin-independent LC wavefunctions (/){x,k'^ = 0.05M^) of pseudo-scalar (solid), scalar 
(short-dashed) and vector (long-dashed) mesons in the chiral limit mo = 0. Mass of the vector meson 
is taken to be a typical value my/2M = 0.9 (see Fig. 4). Rescaled at x = 1/2 just for comparison. 

when the scalar meson appears as a bound state. Then, the difference of the mass between = 0, 
and = 2M greatly affects the shape. Both have peaks at x = 1/2, but widths are significantly 
different. The wavefunction of the scalar meson has a very sharp peak at x = 1/2 (and even divergent 
for k± = 0), which clearly shows the constituent picture. In Fig. 3, we plot (j){x,kj_ = 0.05M^)'s of 
scalar and pseudo-scalar mesons in the chiral limit, which are compared with that of the vector meson 
to be derived soon. Transverse momentum has been taken nonzero in order to avoid the singular 
behavior at k± = 0. The wavefunctions are rescaled at x = 1/2 just for comparison. 

5.3 Vector mesons 

We have been concerned if the bound-state equations for the transverse and longitudinal modes 
(Eqs. (4.24) and (4.26)) are equivalent to each other. This can be explicitly shown as follows. First 
of all, we can verify that both the equations (4.24) and (4.26) derive the same equation for a vector 
meson mass my with the extended parity invariant cutoff (5.4). Namely, Eq. (5.5) for the transverse 
and longitudinal potentials 



leads to the same equation. This, together with the fact that the spin-independent LC wavefunctions 
for transverse and longitudinal modes have the same functional form (5.1), ensures the equivalence of 
the two bound-state equations. If the masses of transverse and longitudinal vector mesons coincide with 
each other, then so do the spin- independent LC wavefunctions, as well as the bound-state equations^'^ . 

^■^Notice that we can derive the eigenvalue equations simply by inserting the LC wavefunction (5.1) into the (r.h.s. 
of the) bound-state equations with the extended parity invariant cutoff. Since the LC wavefunction (5.1) is a direct 
consequence of the bound-state equations, to obtain the same equation for mv means that the original equations are 
also equivalent to each other. 



gy pI + M^ -2y{l-y)pl 
(27r)3 y2(i_y)2 



(2^)3 y{l-y) 
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Indeed, after some calculation, one can see that Eq. (5.5) with i = T and i = L reduces to the same 
equation for the vector meson mass my. 
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(5.15) 



where we have defined a dimensionless coupling constant gy = (A^/47r^) gy and ry similarly as for 
the scalar and pseudo-scalar cases: 

^v-(S)^ (5.16) 

This completes the proof of the equivalence between the transverse and longitudinal equations. 

Let us consider the physics consequences of the eigenvalue equation (5.15). A physical bound-state 
should appear only when the ratio ry is in the range < ry < 1. Equation (5.15) has a solution in 
this region when the strength of the coupling constant ^y is in the range ^y"™'* < 5y < 5y°'^^'* defined 

by 



~(min) _ 3 r 2m f 1 + 



~(m£ix) 

9y 



1 



(5.17) 



Two limiting cases = g^ and gy = gy correspond to ry = 1 (loose binding limit), and ry = 
(tight binding limit), respectively. When M/A ^ (/? ^ 1), the physical bound-state region shrinks 
~(min) _^ y^"^^^ ^ while it becomes wider as M/A grows large. The existence of 5y^™^ is consistent with 
the observation that there is no bound state in the NJL model without the vector interaction. Similar 
behaviors have been found in Ref. [19]. 

One can find approximate solutions to Eq. (5.15) when the vector coupling gy is close to either 

. . —(mill) . —(max) , 

its minimum gy or maximum gy . iNameiy, 



ry ~ 



[ (g^"^)-^ - (gy)-^ 
1 7r(l-/32) 

(^v)-^ - (6^S"^V^ 



(1-^2) (lni±f-2^) 



P ~ ^ —(mill) 

for gy <,gy ' , 

c ~ ^ ~(max) 

for gy ^ gy ■ 



(5.18) 
(5.19) 



When we derived the bound-state equations (4.24) and (4.26), we picked up the leading contribution 
in the (integral form of the) eigen- value equation for my (see Appendix F). However, the solution to 
the resulting eigen-value equation (5.15) has highly nontrivial dependence upon gy. This is evident 
in the above approximate solutions to the equation. 

In Fig. 4, a numerical solution to Eq. (5.15) is shown as a function of gy, where the constituent 
quark mass is taken to be M/A = 0.4 {(3 = 0.82) as an example. As we expect, a bound state appears 
for gy larger than the critical value ^y""^"*, and as one increases gy, the mass starts to decrease from 
the threshold value 2M. The value of critical coupling constants are exactly the same as the values 
predicted by the analytic calculation. When (5 = 0.82, they are gy^™^ = 0.95 and 5y"^^^ = 2.72, which 
are indicated as two vertical dotted lines in the figure. Besides, we have plotted the approximate 
solutions Eqs. (5.18) and (5.19) in the same figure, and they are in good agreement in the regimes 
where each approximation is valid. 

It is also interesting to see my as a function of gy for different values of the constituent mass. 
This is shown in Fig. 5. Unlike in Fig. 4, the vertical axis is my/A instead of my/2M since we 
change the value of constituent mass M/A = 0.1 ~ 0.5. Though the eigenvalue equation (5.15) does 
not explicitly depend upon and mo, they indirectly affect the vector meson's mass through the 
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Figure 4: Vector meson's mass mv/2M as a function of ^y. The numerical solution to Eq. (5.15) 
(solid line) is compared with its approximate solutions Eq. (5.18) (short-dashed) and Eq. (5.19) (long- 
dashed). Constituent quark mass is taken to be M/A = 0.4 (see Fig. 1). A bound state appears only 
in the regime ^y"™^ = 0.95 < gy < g^^^^^ = 2.72, which are indicated by two vertical dotted lines. 
The horizontal dotted line corresponds to the threshold. 



my/ A 




Figure 5: Vector meson's mass my /A as a function of ^y for different values of constituent mass M/A. 
From bottom to top, M/A = 0.1, 0.2, 0.3, 0.4, 0.5. The line for M/A = 0.4 is the same result as shown 
in Fig. 4. Below the dashed line is the physical region. These results can be re-interpreted as either 
'^o or dependences from Fig. 1, since the mo and dependences enter the vector meson's mass 
only through the constituent mass M. 
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constituent mass M (cf. Fig. 1). Thus, to study the M dependence of my also imphes to study the 
or niQ dependence. The result in Fig. 5 is consistent with the 3 dimensional plot of my/ A as a 
function of and in the chiral limit, which was presented in the letter [10]. As M/A is increased 
(i.e., going to larger ^g), the bound-state region becomes larger as was discussed around Eq. (5.17). 

Lastly, we compare the spin-independent LC wavefunction of the vector meson with those of 
the scalar and pseudo-scalar mesons. This is shown in Fig. 3 as already mentioned. We plot three 
wavefunctions for nonzero transverse momentum k'^ = 0.05M^ in the chiral limit mo = 0. Mass of 
the vector meson is taken to be a typical value my /2M = 0.9 (see Fig. 4). All the wavefunctions 
are normalized at x = 1/2 for comparison. As we discussed before, the difference of shape is due to 
different values of the mass. In the chiral limit, we have a relation = m-,^ < uiy < ma = 2M. Thus it 
is reasonable that the wavefunction of the vector meson lies in between those of the scalar and pseudo- 
scalar mesons. The constituent picture works better for a vector meson than for a pseudo-scalar 
meson. 

6 Summary and discussions 

We have presented a detailed analysis of the application of light-front quantization to light mesons. 
Our framework is one of the most straightforward approaches towards relativistic bound-state systems, 
and indeed allows us to obtain the light-cone wavefunction and mass of the meson simultaneously by 
solving the Hamiltonian eigenvalue equation. One of the subtleties in treating light mesons in the 
light-front formalism is the chiral symmetry breaking. However, we know how to describe it at least 
in the Nambu-Jona-Lasinio model as discussed in the present paper. All the calculation was done 
within the NJL model with the vector interaction and in the large A'^ limit, where a meson is described 
as a composite state of a quark and an antiquark. We mainly focused on the problems of describing 
vector mesons and understanding their properties, though we also reproduced the previous results for 
the scalar and pseudo-scalar mesons in Rcf. [4]. The vector mesons can be formed in the NJL model 
if one adds the vector interaction which strengthens the attractive force between the quark and the 
antiquark. For the light mesons, the effect of chiral symmetry breaking is of special importance. We 
followed the same procedures as in Ref. [4] to incorporate this effect. We are able to obtain the LC 
Hamiltonian in the broken phase, which then can be used in deriving the bound-state equations. In 
the NJL model, we have to be careful about the regularization scheme, but once it is properly done, 
we can obtain reasonable results about the LC wavefunctions and masses of light mesons. Mass of 
the vector meson decreases as one increases the strength of the vector interaction. This behavior is 
consistent with Refs.[19, 20] which also treated the vector meson within the NJL model with the vector 
interaction. 

We utilized several approximations in deriving the bound-state equations. For the vector meson, 
we estimated only the leading effect with respect to the vector coupling. In fact, this was enough 
to form a bound state in the vector channel, but somewhat unexpected behavior of the mass that 
it vanishes at some critical vector coupling ^y^^'^^ is probably due to this approximation. For the 
scalar and pseudo-scalar mesons, we ignored the effects of vector interaction, because bound states 
are possible even with the scalar interaction. We know that the vector interaction can modify the 
structure of pseudo-scalar mesons, and it would be important to study this effect. 

The methods we developed in the present paper seems specific to the NJL model. This is not 
true. Actually, the applicability of the formalism is not limited to the NJL model, but is much wider. 
For example, the same formalism can be applied to any fermionic models without gauge fields, which 
include models with nonlocal current-current interaction generated by either one gluon exchange or 
instantons. 

There remain several interesting topics to be considered in the future. First, now that we obtain 
the LC wavefunctions of light mesons, we can use them to compute the physical observables such as 
the distribution amplitudes. Of special importance is the physical form factors. Such analysis was 
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indeed done for a pseudo-scalar meson [15], but a similar analysis is possible for the vector meson [21]. 
Next, our framework can in principle treat the axial vector state in a similar way, and it would be 
interesting to investigate its bound-state structure. 
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A Notations 

Our notation is the following: The light-cone coordinates are = {x^ ±x^)/\/^, and x\_ = = 
x±, which are written together as x^^ = {x~^ , x~ , x\_) . The derivatives are d± = d/dx^ and d± = 
d/dx±. We define the spatial vector x for coordinates by x = {x~,x\^. On the other hand, for 
momenta, wedefinep = {p'^ ,p]_,p'^) (notice that p"*" = p- is conjugate to x~) and pa; = —p'^x~ -\-p\x^j_ 
so that px = p^x^ + p^x^ — pj_Xj_ = p^x^ — px. We use /U, v for Lorentz indices of four vectors, i,j 
for transverse coordinates 1,2, and a,P for spinor indices. 
The projection operator A± is defined by 

A±^^(1±7V) (A.l) 
and satisfies (A±)^ = A±, and A+-|-A_ = 1. Some useful formulas are A_7+ = 7"'"A+, A+7~ = 7~A_. 



B LF Quantization of a free fermion: spinor basis 

Here we summarize the LF quantization of a free fermion based on the plane wave solutions to the 
Dirac equation. This gives our "spinor basis" which is useful when we treat the case with vector 
interaction. 

Consider a classical plane wave solution = u{p) e^^^^ to the free Dirac equation: 

(7>^-M)u(p) = 0, (// = +, -,1,2) (B.l) 

where the four momentum p^^ must satisfy the dispersion relation 

So far, the LC energy p" can take either positive or negative value, but, we redefine the negative energy 

solution u{—p^, —p~^, —pi_) {p~ > 0, p"*" > 0) by v{p~ ^p'^ ^p_[_), and we always treat p~ to be positive. 
We call u{p) and v{p) the positive and negative energy solutions, respectively. It should be noticed 
that the longitudinal momentum p+ is also positive, because the sign of the LC energy and that of 
the longitudinal momentum is correlated as is shown in Eq. (B.2). The "positive energy solution" 
u{p~^ ,p±, s) satisfies = {^^Pfj, — M) u{p~^,pj_, s) and u{p~^ ,p±, s) {j^Pn — M) = 0, where s is an index 
to distinguish two different solutions. Similarly, the "negative energy solution" v{p'^ ,p±, s) satisfies 
= {—'^^Pfi — M) v(j)'^ ,p±, s) and v(j)'^ ,p±, s) (— 7'^P/i — M) = 0. Normalization of the solutions are 
given by 

u{p, s) u{p, s') = -v{p, s) v{p, s') = 2M 5ss' , 
and the foUowings are useful: 

u{p,s)j'^u{p,s') = v{p,s)j'^v{p,s') = 2p^'Sss', 
Y^u{p,s)u{p,s) = p^-^i^ + M, 

s 

^v{p,s)v{p,s) = p^Y-M- 
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The fermion field can be expanded by use of these complete solutions, 

2 roo POO 

where b(p^,p±, s) and S{p'^ ,p±., s) arc c-number coefficients. Defining ■0 = A^ip + A^tp = + i/j^, 
the good component ^+ is expanded as 



ip+{x) = ^ / dp^ / d'^p± — 



E [b{p+,p^,s) u+ip+,p^,s)e-'''-^ + d\p+,p^,s) v+ip+,p^,s)e'P--] , (B.4) 



X 

s=± 



where u+ = A^u and u+ = A+v. These projected spinors satisfy the following relations: 

ul {p, s')u+ (p, s) = v^ip, s')v+ {p, s) = V2p-^ dss' 

Yl (p^ ^)'"+ (p^ ^) = 12'"+ ^p^ ^p^ = ^^H-p"^ ■ 

s s 

Since ip^ is a dependent variable, we impose the quantization condition only on the good compo- 
nent ^+ (see Eqs. (2.5), (2.6)). The mode expansion for ijj^ is given in the same form as Eq. (B.4), 
but now b{p'^ ,p±, s) and d^{p'^,p±,s) are operators. 

The mode expansion of the bad component tp^ should be given through the constraint equation 
il;^{x) = 2i5r(^7l^* + A'/)7+^+ by using Eq. (B.4). One can explicitly show that V' = V'+ + V'- 
reproduces the mode expansion Eq. (B.3). 

C Boson expansion method as the expansion of M{p,q) 

In the boson expansion method, we rewrite the bilocal operators M(p, q) defined by Eq. (3.1) with re- 
spect to bosonic operators so that the commutation relation between M(p, q)'s is correctly reproduced. 
Physically this corresponds to capturing a composite state made of a fermion and an antifermion, as a 
bosonic state with additional "structure" due to statistics. Especially, when N is large, the composite 
state behaves as a boson, and thus large N expansion of the bilocal operator is naturally realized by 
the boson expansion method (of the Holstein-Primakoff type) . 

The bilocal operator M{p,q) satisfies the following commutator: 

: Maia2(Pi,P2) :> : Ma^a2(Qi>92) : ] (C-l) 
= : Ma^p,^{p^,q2) : (A+)„,/3^(^(p2 + Qi) - : Mp^^^{qi,p2) : (A+)^2ai^(Q2 + Pi) 

+ N{A+)g.^ai^{q2 + Pi){A+)^.^0^6{p2 + qi) 

X {e{pi)e{pt)e{-qt)ei-q+) - e{-pt)ei-pt)eiqt)e{q+)), 

and its vacuum expectation value is 

(0|M„^(p, q)\0) = N9{p+)6{p + q)(A+)^,. (C.2) 

This complicated commutator can be reproduced by replacing M (p, q) by some simple functions of a 
(bilocal) bosonic operator. 
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Let us first introduce bilocal bosonic operators Bg^p(p,q) which have Dirac indices. We impose 
they satisfy the following commutation relations: 

[Baia2{Pi,P2), Bl^f3^{qi,q2)] = {A+)a,f3,S{pi - q^){A+)p^aJ{q2 - P2), (C.3) 

[5aia2(Pl,P2),%/32(9l,Q2)] = > 0), (C.4) 

and the following properties 

(A+)ai7i-B^l72(Pl)P2) = -^«i72(Pl)P2)) -^7172 (Pl ) P2) (■^+ )72a2 =-^7102(^1)^2); (C5) 
(A+)7iai^^l72(Pl,P2) = Bi^^^{Pi,P2), St^^^(Pi,P2)(A+)a272 = B\^^^{Pi,P2) , (C.6) 

and so on. Then the operators : M : can be represented as 

poo roo 

: M-+2(Pi'P2) : = / dq+ (fq^Bl^^{-p^,q)B^^^{p2,q) 
Jo J— 00 

= Aa2aiiP2,-Pl), (C.7) 

roo POO 

:M+-,(Pi,P2): = - / dq+ d\±B\^^{q, -p2)B^aAQ,Pi). (CS) 

Jo J —00 

: M++,(Pi,P2) : = / dq+ d^q±WN - A)a2'r{P2,Q)BjaAQ,Pi), (C.9) 

Jo J —00 

roo roo 

:M--a,(Pi,P2): = dq+ d\^B\^^{q,-p2){VN^)^^,{q,-Pi), (CIO) 

Jo J— 00 

where the upper indices stand for the signs of the longitudinal momenta. This is the boson expansion 
method of the Holstein-Primakoff type, represented in the "Dirac basis" (namely, the bosonic operators 
explicitly have the Dirac indices). This is essentially the same as the one used in Ref. [4]. Note that 
the bosonic operator is of 0(iV°), and there is explicit dependence on N on the right-hand sides. Thus, 
one can expand the right-hand sides with respect to 1/^/N, yielding Eq. (3.6) in the text. It should 
be noticed that the natural expansion parameter is instead of 1/N because we have expanded 

^/l — A/N . The first few terms are given as follows: 

m^^lip.q) = {A+)p^6ip + q)e{p+)9i-q+) , (C.ll) 
m(^'j(p,qr) = + ^^^(-p, , (C.12) 

mgJ(p,Q) = l[dk]Bi^{-p,k)Bp^{q,k)e{-p+)9{q+) 

- I [dk] B\^{k, -q)B^a{k,p)9{p+)9{-q+) . (CIS) 

The leading term m^^^ {p, q) comes from the vacuum expectation value of M (p, g) . The result m^^^ (p, q) 
means that the bilocal operator M (p, q) can be treated as a bosonic operator in the first nontrivial 
leading order. 

One can further introduce bilocal bosonic operators B{pi,si : P2, S2) in the "spinor basis", namely, 
with spins instead of simple Dirac indices. The relation between operators in the "Dirac basis" and 
in the "spinor basis" is given by 

Ba^a2{Pl,P2) = —f^— ^ X! "+"l(f'l'*l)^+a2(P2>S2)-B(Pi,Sl :P2,S2), (C.14) 

where the good components of the spinor = A+u and v-\- = A+v are introduced in Appendix B. 
The operator B{pi, si : P2, S2) satisfies a very simple commutator (as shown in the text): 

[-B(Pi, si : P2, S2), B\p[, s[ : P2, 4)] = ^sis' ^(P^ " Pi)^s2s' ^(P^ ~ ^^2) • 
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D Higher order solutions to the fermionic constraint and the LC 
Hamiltonian in the case without vector interaction 



D.l Higher order solutions 

The bilocal fermionic constraints (3.4) and (3.5) for the higher orders can be cast into a compact form 
similar to that of the lowest order: 



where and are known functions expressed by m(") and lower order quantities s*^"*), p^"*) and 
m^"^) with m < n. For example, in the next leading order 0{-s/N) with n = 1: 

and in the next-to-next leading order 0{N^) with n = 2, one finds 

f^^Hp,i) = (^'y-^ {.g(p,,)-.g)(,,p)} 

'^■/^{->^-^-0p^+n5p(^)]^^(fc,0 

(fe,Om(;j(g-fc-Z,p)| , (D.3) 



2q+ 



a/3 



and similarly for G^'^\ 

The integral equation (D.l) has a very simple structure. One can solve it by integrating the 
equation so that the integral on the left-hand side becomes the same as the second term in (D.l). 
Namely, the solution is 

sW(p,q)\ _ /F(")(p,q)\ e(p+) ^ /FW(fc,p + g-fc)\ 

p(") (P, q)J~\ G^") {p, q) ) ^^"^ q+ J-oo (27r)3 {k,p + q-k)J ^ ^"""^^ 

where Kg is a constant defined in Eq. (3.18). Since and are expressed by m^") and lower 
order quantities, we can in principle write down the solutions for any n, order by order. 

D.2 Deriving the LC Hamiltonian 

The LC Hamiltonian is the evolution operator in the LC time (a;+) direction. In deriving the LC 
Hamiltonian, we need to follow Dirac's procedure for constrained systems regarding x'^ as time, since, 
as we already saw, we have the fermionic constraints on the light front. 
Then, one finds the following hermitian Hamiltonian: 



4/ 



mfdi A-^' - {di^)iY A+* - mo*^' , (D.5) 



where i = 1,2 and we have used the fermionic constraint (2.3) since the model is of the second class (in 
Dirac's terminology). If one takes its face value, this Hamiltonian, being independent of the coupling 
constant Gs, looks equivalent to the free Hamiltonian. However, if one rewrites it with respect to 
the physical degree of freedom, - namely the good component of the spinor -, then one obtains a 
non-trivial Hamiltonian that depends explicitly on the coupling constant. This is of course because 
the bad component ijj- is subject to the fermionic constraint (2.3) which carries the information of 
interaction. 
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Since the Hamiltonian (D.5) is bilinear with respect to the spinor, one should be able to rewrite 
it by the bilocal operators. This is indeed the case if one introduces a new bilocal operator V| (see 
Eqs. (3.14) and (3.15)). As mentioned in the text, the bilocal operator Vj can be related to the 
known ones, M, S, P, through an equation which is easily derived from the fermionic constraint. In 
momentum space, it is given by 

iq^y'Mq) = -^{7i(7i9i + "io)}^^M„^(p,q)+ ^{(7igi + mo)7i}^^M„;3(qf,p) 
+ ^ / [MapiP, g - fe - {tI (Sr + i75 Pr) {k, I) 



{(Sr - Z75 Pr)71}^^ {k,l) Mapiq -k-l,p) 



(D.6) 



Therefore, given the solutions Sp*^ and p^i , one can immediately obtain vj the expansion coefficients 
of VJ. Hence, one can compute the LC Hamiltonian order by order. 

(n) (n) 

We substitute the "broken" solutions of s^' and pj^' to obtain the Hamiltonian which describes the 
broken phase of the chiral symmetry. It turns out that the lowest order /tg°^ is divergent but is just 
a constant, and thus can be neglected for the present purpose. The next order ^^\p,—p) and 
s^\p,—p) become strictly zero, and we have /ig^^ = 0. Therefore, the first nontrivial contribution 
is given by the next-to-leading order n = 2. After straightforward calculation, one arrives at the 
following expression: 

= -/.l^P^n,S(p,-p) 

g^Ks f°° (fp(Pq(fk(fl,, , 

X {F("(p,Q)FW(fe,() + 0("(p,g)0("(fc,i)} , 

where we have again ignored a c-number contribution, which is unimportant for our present purpose. 
It is convenient to express the functions F^^^ and G^^^ given in Eq. (D.2) by using the spinors (below, 
p+ > 0, g+ > 0): 

F^^\-p,-q)+F^^\-q,-p) = -j^== ^ ^ {u{p,si)v{q,S2)}B^{p,Si : q,S2) + (p ^ q), 

G(i)(-p,-q) + G(^)(-g,-p) = _L_ J2 ^ {u{p,si){ij5)v{q,S2)}B\p,si : q, S2) + (p ^ q) . 

V Q si,s2 v2 

Using these expression, together with the explicit representation for m(^'^)(p, q) (cf. Eqs. (C.12), 

(2) 

(C.13)), one can rewrite hg with respect to the bilocal boson operators to obtain the final expression 
(3.17). 



E Meson states with nonzero transverse momentum P± 

Here we present the case with nonzero transverse momentum Pj_ 7^ which requires only a small 
generalization of the case with P± = 0. In fact, this is almost trivial from the viewpoint of the light- 
front formalism and even not necessary for the present purpose of this paper. However, it is useful 
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to have a formalism with nonzero transverse momentum when we wih compute the form factors of 
mesons [21]. 

A generic meson state with nonzero transverse momentum 7^ is given by 



meson 



■P) = P+ dxrdx2 5{l-Y,Xi) (fki±cfk2±6{J2ki±) 
Jo ^-00 .^^^ 

Y^^{xi,ki±,si,S2)B\p-^,si:p2,S2) |0) , (E.l) 



X 

S1,S2 



where Pi and P2 are the momenta of a quark and an antiquark P = p^ + P2, and the relations to the 

"relative" momentum coordinates Xi and ki± are given by p^ = (xiP~^ , XiP± + kij_). These variables 
Xi and ki± are invariant under any boost transformations. If one performs X2 and k2± integrations, 
one obtains 

rl roc 

I meson ;P) = P+ / dxi d^ki± ^ ^xi,ki±, si, S2)B\p^, si : P - p^, S2) \0) , (E.2) 

"'O J-°° S1,S2 

As was done in the case with P± = 0, the kinematical structure of the LC wavefunction ^{xi,ki±,Si) 
is determined through the interpolating field of each meson. For example, the LC wavefunctions of 
the scalar and pseudo-scalar mesons are written as follows: 

$^(xi,fci_L,Si) = (j)nixi,ki±) ^ =(u{p^,Si)ij5v{P -pi,S2)\ , 

^aixi,ki±,Si) = (f)aixi,ki±) ^ = (u{pi,Si)v{P - Pi,S2)\ . 

^\/Xl[l Xl) 

In fact, one can explicitly show that the right-hand sides are independent of the total momentum P. 
Therefore, the resulting bound-state equations for the spin-independent part of the LC wavefunctions 
are exactly the same as those derived in the frame with P± = 0. 

What is less trivial is the case of the vector meson due to the projection onto the three physical 
modes. From the interpolating field, one obtains the expected result: 

^pixi,ki±,Si) = <^p(xi,fci±) ^ e^(A,P){u(pi,si)7^t;(P-pi,S2)} , (E.3) 

^\/Xl[L Xi) 

but the polarization vector in a generic frame becomes more complicated [13] 

e''(A = ±l,P) = (0, 



:fP^-iPy Tl 



V2P+ 'V2'V2'' 

e>'{X = 0,P) = (— 

my 2 my P^ my my 

satisfying e^(A, P)e'*(A', P)* = —Sxy, P^e^^{\,P) = 0. Of course this vector reproduces Eq. (4.13) 
when Px = Py = 0. It is very important to recognize that the spin-dependent part of the LC 
wavefunction, eju(A, P) u{pi, si) v{p2, S2) with p^ = {xiP~^,xiP± + ki±) and P2 = {x2P~^,X2P± + 
^2±)) does not depend on the total momentum P. This can be explicitly verified as follows. Using 
the Kogut-Super convention, one finds for A = -|-1 

€^( + 1, P) W((X1P+, XlPx + fcix), .Si) r viix2P+,X2P± + k2±), S2) 

V2 ( -M xi{kl^ + ikl^)\ 

./x^\x2{kl^ + ikl^) j^^^^ ' 

and similarly for other polarizations. Therefore, from these explicit calculations, one concludes that 
the bound-state equations (for the relative motion of a quark and an antiquark) in the frame with 
nonzero transverse momentum are the same as those with Pi =0. 
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F Bound-state equation of longitudinally polarized vector meson 

In this Appendix, we first derive the potential (4.25) for the longitudinally polarized vector meson 
and then the bound-state equation (4.26) after taking the leading contribution with respect to the 
coupling constant ^y. 

The bound-state equation for the longitudinal vector meson is derived from Eq. (4.20) where the 
Hamiltonian and the state are given by Eq. (3.27) and Eqs. (4.3), (4.7), (4.15), respectively. Explicitly, 
for the longitudinal mode (we use rrii^ for the vector meson's mass), 



mi 

2P+ 1 2P+ ' ' mL ' j 1 2P+ rriL 



where we have chosen the frame P = {P^,P±) = {P^,0±). The following relation is useful when we 
evaluate this equation: 

1 k'^ + 

u{k, si)7~t;(P - k, S2) = ~(^Tj2 2a^(l - x) ^^^^'^'"^^ ~ ^2)- 
Then, a straightforward calculation yields the following bound-state equation: 

I TTT-L X(l — X) J 

= |' + <^(r^|i^(r^^L^^'^^) 

2gy kl + M^ /•! 



(27r)3 x(] 



dy j dH^ {{ll + m2) + mly{l - y)} 



(27r)3yo "J ^^'^ ' miy{\-y) y{\ - y) 

- (r -7 ^^'^ ^) r / ^^'^ ^) ^ ^ } • 

Note that this result is independent of the total longitudinal momentum P^ . It is not difficult to 
identify the origin of each term. On the right-hand side, the first term comes from the kinetic term of 
the Hamiltonian, the second and third terms are from and the last term is from uyi. 

Remarkably, this complicated expression can be greatly simplified if one recognizes that the equa- 
tion is made of three different types of integral. These constants are defined by the following: 

C ^ j\yJd'l^{mly{l-y) + {ll + M')}^^, 

V = dy d b± 2 n V ' 

Jo J mly{l-y) 

In fact, these three constants are not independent of each other. One can verify this by first multiplying 
Eq. (F.l) by x{l — x) and then integrating it over x and k±: 

u = 5 % + 2-%rir]. (F.2) 

l-2j^j^dxjd^k^ (27r)3^' 
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This relation helps to simplify the third term on the right-hand side of Eq. (F.l). Indeed, substituting 
this into the third term, one gets 



x{l — x) 



= 1-2 



(27r)3 Jo 



dz J (fq± 



9y 



4 {kl + M2 



(27r)3^ (27r)3 ^ ml x{l - x) - {kl + M^) \ (2 
Plugging the definition of ^ back to this equation, one finally obtains 



(fk± 



"'L + ^h — 7 

x(l — x) 



(27r)3 



2 + 



6l(x,A;_l) 

4(A:2 +m2) 



m2x(l-x)-(A:2 +M2) I (27r)3 Jo 



fv 



dx d k± 



This gives the potential (4.25) (notice the plus sign of the second term on the left-hand side while the 
potential is defined with minus sign, see (4.19)). This is still complicated and we will approximate 
this equation by carefully taking the leading contribution with respect to gy 

Let us integrate the above equation over the external variables x, k±. Then, we obtain a simple 
integral equation: 



gy 
(27r)3 Jo 



dx / d^k± 



4 {kl + M2 



m£x(l - x) - {kl + M"^) 



gy 
(27r) 



./;d./.v}. 



This is actually an eigenvalue equation for ttt-l, and we shall treat this equation by taking the lowest 
order contribution with respect to gy. Namely, we treat 



1 = 



gy 
(27r)3 Jo 



dx / d^fc_L 



4 (yt2 + m2) 



ix{l-x) - {kl + M^)' 



m 



In the bound-state equation (F.3), this approximation corresponds to ignoring only the g^ dependent 
term in the square brackets. Then, one obtains 



2 + m2 



<pL(.x,k±) 
• 1 



x{l — x) 



ll + M^ 

y{^-y) 



b^{y,l±). (F.4) 



Integrating this equation again over x, k±, one finds 



dxd k± < mL 



kl + M^' 
x{l — x) 



I gy 
(27r)3 



dxd k±] dydl± 



ll + M^ 



(l)L{yJ±)- 



Using this relation on the right-hand side of Eq. (F.4) finally yields the bound-state equation (4.26). 



2 kj + 



x{l — x) 



Xx,k±) 



gy 
(27r)^ 



dy / 



4(/i + 

y{^-y) 



(f>L(.yJ±)- 
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